Deuteron and antideuteron production in Au+Au collisions at y^s^ = 200 GeV 
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The production of deuterons and antideuterons in the transverse momentum range 
1.1 < pt < 4.3 GeV/c at mid-rapidity in An + Au collisions at y/SNN = 200 GeV has been 
studied by the PHENIX experiment at RHIC. A coalescence analysis comparing the deuteron and 
antideuteron spectra with those of protons and antiprotons, has been performed. The coalescence 
probability is equal for both deuterons and antideuterons and increases as a function of pt, which 
is consistent with an expanding collision zone. Comparing (anti)proton yields : p/p = 0.73 ±0.01, 
with (anti)deuteron yields: d/d = 0.47 ± 0.03, we estimate that n/n = 0.64 ± 0.04. 

PACS numbers: 25.75.Dw 

Keywords: relativistic, heavy ion, collisions, deuteron 

Ultrarelativistic heavy ion collisions are used to study temperature and density, similar to those that existed 
the behavior of nuclear matter at extreme conditions of in the universe a few microseconds after the Big Bang. 



Previous measurements indicate that high particle mul- 
tiphcities [l| and large p/p ratios prevail at the Relativis- 
tic Heavy-Ion Collider (RHIC), which is expected for a 
nearly net baryon free region [2|. As the hot, dense sys- 
tem of particles cools, it expands and the mean free path 
increases until the particles cease interacting ("freeze- 
out"). At this point, light nuclei like deuterons and an- 
tideuterons {d and d) can be formed, with a probability 
proportional to the product of the phase space densities 
of their constituent nucleons |3, Ijl ■ Thus, invariant yield 
of deuterons, compared to the protons [5|, |a] from which 
they coalesce, provides information about the size of the 
emitting system and its space-time evolution. 

PHENIX |3| at RHIC, is a versatile detector designed 
to study the production of leptons, photons, and hadrons 
over a wide momentum range. In this Letter, results on d 
and d production in Au-|-Au interactions at ^/sNN — 200 
GeV are presented. For the sake of brevity, in the rest of 
this Letter, our statements will generally apply to both 
particles and antiparticles. 

The East central tracking spectrometer in the 
PHENIX detector H, Q, "^ is used in this analysis. The 
information from the PHENIX Beam-Beam Counters 
(BBC) and Zero-Degree Calorimeters (ZDC) is used for 
triggering and event selection. The BBCs are Cerenkov- 
counters surrounding the beam pipe in the pseudorapid- 
ity interval 3.0 < \r]\ < 3.9, and provide the start tim- 
ing signal. The ZDCs are hadronic calorimeters 18 m 
downstream of the interaction region and detect specta- 
tor neutrons in a narrow forward cone. Particle identifi- 
cation in the central rapidity region is achieved by mea- 
suring momentum (by drift chamber) and time of flight 
(by time-of-flight detector). The drift chamber (DC) and 
two layers of pad chambers (PC) are used for tracking and 
momentum reconstruction [g]. The time-of-flight detec- 
tor (TOF) spans the pseudorapidity range \r]\ < 0.35 and 
A(j) = 7r/4 azimuthally. The TOF consists of plastic scin- 
tillators, with a combined time resolution of « 115 ps. 
The TOF thus provides identification of d and d in the 
transverse momentum (pt) range 1.1 < px < 4.3 GeV/c. 
For pt < 1.1 GeV/c, the signal to background ratio suf- 
fers due to multiple scattering and energy loss effects. 

The dataset for this analysis includes 21.6 million min- 
imum bias events. The minimum bias cross section cor- 
responds to 92.2^3^% of the total inelastic Au-f Au cross 
section (6.9 b) |9j. Using the momentum determined by 
the DC, which has a resolution of Sp/p ~ 0.7% © l%p 
GeV/c, and the time of flight from the event vertex pro- 
vided by the TOF, the mass of the particle is deter- 
mined. The d and d yields are obtained by fitting the 
mass squared distributions to the sum of a Gaussian sig- 
nal and an exponential background. Examples of mass 
squared distributions with fits for antideuterons in mini- 
mum bias collisions are shown in Fig. ^ 

The raw yields are corrected for effects of detector 
acceptance, reconstruction efficiency and detector occu- 
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FIG. 1: (color online) Histograms of the mass squared for 
identified antideuterons in the transverse momentum range 
1.1 <pt < 3.5 GeV/c (in 400 MeV/c increments), with Gaus- 
sian fits including an exponential background. 



pancy. Corrections are determined by reconstructing sin- 
gle deuterons simulated using GEANT [lOj and a detec- 
tor response model of PHENIX, using the method de- 
scribed in '6J. The track reconstruction efficiency de- 
creases in high multiplicity events because of high de- 
tector occupancy. This effect can be slightly larger for 
slower, heavier particles, due to detector dead times be- 
tween successive hits. Occupancy effects on reconstruc- 
tion efficiency (« 83.5% for 0-20% most central events) 
are evaluated by embedding simulated single particle 
Monte Carlo events in real events. Since the hadronic 
interactions of nuclei are not treated by GEANT, a cor- 
rection needs to be applied to account for the hadronic 
absorption of d and d (including annihilation). The d- 
and d-nucleus cross sections are calculated from parame- 
terizations of the nucleon and anti-nucleon cross sections: 



'^d/d,A — [^/^7iJN^+ ^d] 



(1) 



The limited data available on deuteron induced in- 
teractions lllj indicate that the term A^ is indepen- 
dent of the nuclear mass number A and that A^ = 
3.51 ± 0.25 mb^/^. The hadronic absorption varies only 
slightly over the applicable pr range and is « 10% for 
d and ~ 15% for d. The background contribution from 
deuterons knocked out due to the interaction of the pro- 
duced particles with the beam pipe is estimated using 
simulations and found to be negligible in the momentum 
range of our measurement. 

Figure 121 shows the corrected d and d invariant yields 
as a function of transverse mass {iriT in the range 1.1 < 
Pt < 4.3 GeV/c, for minimum bias events, and two cen- 
trality bins: 0-20% (most central), 20-92% (non-central). 
The 20-92% centrality bin is dominated by mid-central 
events, due to larger track multiplicities relative to pe- 
ripheral events. 
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FIG. 2: (color online) Corrected spectra for deuterons (left 
panel) and anti-deuterons (right panel) for different centrali- 
ties are plotted vs rriT- Error bars indicate statistical errors 
and grey bands the systematic errors. Values are plotted at 
the "true" mean value of mr of each bin, the extent of which 
is indicated by the width of the grey bars along x-axis. 



TABLE I: The inverse slope parameter T^ff obtained from 
a niT exponential fit to the spectra along with multiplicity 
dN/dy and mean transverse momentum (pr) obtained from 
a Boltzman distribution for different centralities: 



Teff [MeV] 


Deuterons 


Anti-deuterons 


Minimum Bias 


519 ± 27 


512 ± 32 


0-20% 


536 ± 32 


562 ± 51 


20-92% 


475 ± 29 


456 ± 35 


dN/dy 



Minimum Bias 

0-20% 

20-92% 



0.0250 ± 
0.0727 ± 
0.0133 ± 



0.0006(3tat.) 

0.005(3ys.) 

0.0022{stat.) 

Q.0141(sys.) 

0.0004(stat.) 

0.0029(sys.) 



0.0117 ±: 

0.0336 ± 
0.0066 ± 



0.0003(stat.) 

0.002(sj/s.) 

0.0013{stat.) 

Q.0057(3ys.) 

0.0002(stat.) 

0.0015(sK3.) 



{pt) [GeV/c] 



Minimum Bias 

0-20% 

20-92% 



1.54 ± 
1.58 ± 
1.45 ± 



0.04,(atat.) 
O.ialays.) 
0.05{stat.) 
O.ialays.) 
0.05(stat.) 
0.15(sys.) 



1.52 ± 
1.62 ± 
1.41 ± 



O.Q5{stat.) 

0.12(aya.) 

0.07{stat.) 

0.1(sys.) 

0.06(stat.) 

0.15(sys.) 



Systematic uncertainties have several sources: errors 
in particle identification, DC-TOF hit match efficiency, 
the uncertainty in momentum scale, d and d hadronic in- 
teraction correction, and uncertainty in occupancy cor- 
rections. All the systematic uncertainties are added in 
quadrature, depicted by grey bars in Fig. [3 

The PT spectra Ed^N/d^p are fitted in the range 
1.1 < Pt < 3.5 GeV/c to an exponential distribution in 
TTiT — \Jpt + "^ ■ The inverse slopes (T^ff) of the spec- 
tra are tabulated in Tabled The deuteron inverse slopes 
of Tf.ff = 500-520 MeV are considerably higher than the 
Te// = 300-350 MeV observed for protons [iQ. The in- 
variant yields and the average transverse momenta {{pt}) 
are obtained by summing the data over pT and using 



a Boltzmann distribution: 



d^N 



lizniTdrnTdy 



oc rriTe 



-rriT/T^jf 



(which gives a slightly better -x^ /n.d.f. — 4.8/3 vs. 
X^ /n.d.f. = 5.6/3 for exponential fit) to extrapolate to 
low rriT regions where we have no data. The extrapo- 
lated yields constitute w 42% of our total yields. The 
rapidity distributions, dN/dy, and the mean transverse 
momenta, (pt), are compiled in Table Q] for three differ- 
ent centrality bins. Systematic uncertainties on dN/dy 
and (pt) are estimated by using an exponential in pT 
and a "truncated" Boltzman distribution (assumed flat 
ior Pt < 1-1 GeV/c) for alternative extrapolations. 

With a binding energy of 2.24 MeV, the deuteron is 
a very loosely bound state. Thus, it is formed only at 
a later stage in the collision, most likely after elastic 
hadronic interactions have ceased; the proton and neu- 
tron must be close in space and tightly correlated in ve- 
locity to coalesce. As a result, d and d yields are a sensi- 
tive measure of correlations in phase space and can pro- 
vide information about the space-time evolution of the 
system. If deuterons are formed by coalescence of pro- 
tons and neutrons, the invariant deuteron yield can be 
related [l2| to the primordial nucleon yields by: 



Ed 



d^Nd 



d^Pd 



Bo 



Pd='2pp 






(2) 



where B2 is the coalescence parameter, with the subscript 
implying that two nucleons are involved in the coales- 
cence. The above equation includes an implicit assump- 
tion that the ratio of neutrons to protons is unity. The 
proton and antiproton spectra jQ] are corrected for feed- 
down from A and A decays by using a MC simulation 
tuned to reproduce the particle ratios: (A/p and A/p) 
measured by PHENIX at 130 GeV ^. 
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FIG. 3: (color online) Coalescence parameter B2 vs pr for 
deuterons (left panel) and anti-deuterons (right panel). Grey 
bands indicate the systematic errors. Values are plotted at 
the "true" mean value of pt of each bin, the extent of which 
is indicated by the width of the grey bars along x-axis. 

Figure displays the coalescence parameter B2 as 
a function of pT for different centralities. The de- 
creased B2 in more central collisions implies that in larger 
sources, the average relative separation between nucleons 
increases, thus decreasing the probability of formation 



of deuterons. We also observe that B2 increases with 
Pt- This is consistent with an expanding source because 
position-momentum correlations lead to a higher coales- 
cence probability at larger px- The p-r-dependence of B2 
can also provide information about the density profile of 
the source as well as the expansion velocity distribution. 
It has been shown [1J| that generally a Gaussian source 
density profile leads to a constant B2 with pT as it gives 
greater weight to the center of the system, where radial 
expansion is weakest. This is not supported by our data, 
which shows a rise in B2 with pT- 

Thermodynamic models |4| predict that B2 scales with 
the inverse of the effective volume Ve// (-B2 oc l/Vg//). 
The d and d spectra are affected by radial flow, which 
concentrates the coalescing protons and neutrons, affect- 
ing phase space correlations, thereby limiting the appli- 
cability of a simple thermodynamical model to determine 
an effective source size. 



mined from the particle/anti-particle ratios \l< 
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FIG. 4: (color online) Comparison of the coalescence param- 
eter for deuterons and anti-deuterons {pr = 1.3 GeV/c) with 
other experiments at different values of ^/s. 



Figure ^compares i?? for most central collisions to re- 
sults at lower \/i jH [H [13, El El, 113 • Note that B2 
is nearly independent of -y/s, indicating that the source 
volume does not change appreciably with center-of-mass 
energy (with the caveat that B2 varies as a function of 
Pt, centrality and rapidity). Similar behavior is seen 
for B2 for deuterons [ijj as a function of ^/s. This ob- 
servation is consistent with what has been observed in 
Bose-Einstein correlation Hanbury-Brown Twiss (HBT) 
analysis at RHIC |2l| for identified particles. The coa- 
lescence parameter B2 for d and d, is equal within errors, 
indicating that nucleons and antinucleons have the same 
temperature, flow and freeze-out density distributions. 

The ratio n/n can be estimated from the data based 
on the thermal chemical model. Assuming thermal and 
chemical equilibrium, the chemical fugacities are deter- 



FIG. 5: d/d ratio vs. pt for minimum bias data. The dashed 
lines represent the square of the measured p/p ratio as a func- 
tion of Pt within uncertainties. 

Figure [S] shows that the d/d ratio is independent of 
centrality, and px within errors. The average value of 
d/d is 0.47 ± 0.03, consistent with the square of the ra- 
tio p/p = 0.73 ±0.01 p within statistical and systematic 
uncertainties. This is expected if deuterons are formed 
by coalescence of comoving nucleons and p/p = n/n. Us- 
ing the ratio p/p, the extracted proton fugacity is Ap — 
exp(^p/r) — 1.17 ± 0.01. Similarly, using the d/d ratio, 
the extracted deuteron fugacity is Ad = exp[(/ip + /i„)/r] 
= 1.46 ± 0.05. From this, the neutron fugacity can be 
estimated to be A„ = exp(/i„/r) = 1.25 ± 0.04, which 
results in n/n = 0.64 ± 0.04. These estimates, along with 
equality of B2 for d and d indicate that, within errors, 
Mn > Mp- This is expected since the entrance Au-|-Au 
channel has larger net neutron density than net proton 
density. 

To summarize, the transverse momentum spectra of d 
and d in the range 1.1 < px < 4.3 GeV/c, have been 
measured at mid-rapidty in Au+Au collisions at ^/snn 
— 200 GeV, and are found to be less steeply falling than 
proton (and antiproton) spectra. This behavior is con- 
sistent with a constant (flat) source density profile. The 
extracted coalescence parameter B2 increases with px, 
which is indicative of an expanding source. B2 decreases 
for more central collisions, consistent with an increasing 
source size with centrality. The B2 measured in nucleus- 
nucleus collisions is independent of -y/sjvw above 12 GeV, 
consistent with Bose-Einstein correlation measurements 
of the radii of the source. B2 is equal within errors for 
both deuterons and anti-deuterons. From the measure- 
ments, it is estimated that n/n = 0.64 ± 0.04. 
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